Calculation of parity nonconservation in neutral ytterbium 
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We use configuration interaction and many-body perturbation theory techniques to calculate spin- 
independent and spin-dependent parts of the parity nonconserving amplitudes of the transitions 
between the 6s 2 1 So ground state and the 6s5d 3 Di excited state of 171 Yb and 173 Yb. The results 
are presented in a form convenient for extracting spin-dependent interaction constants (such as, e.g., 
anapole moment) from the measurements. 

PACS numbers: 11.30.Er; 31.15.A- 



I. INTRODUCTION 



The use of atomic ytterbium to study parity noncon- 
servation (PNC) in atoms was first suggested by De- 
Mille [l|. The measurements are in progress at Berke- 
ley and afters years of hard work [2|-|5| the first results 
of PNC measurements are finally reported [|| . As it was 
expected the PNC in ytterbium is strongly enhanced, 
being two orders of magnitude larger than in cesium Q. 
The cesium PNC experiment together with its interpre- 
tation [sMTolj in terms of nuclear weak charge provides 
the best current atomic test of the standard model (see, 
also a review [ll| ) . It is also the only measurement of the 
nuclear anapole moment which is produced by the PNC 
nuclear forces }12| . The extraction of the weak nuclear 
charge from the PNC measurements relies on atomic cal- 
culations. The interpretation of the PNC measurements 
in ytterbium similar to what was done for cesium is not 
possible due to limitations of atomic theory. Ytterbium 
has complicated electron structure and calculations for it 
on the same level of accuracy as for cesium are not pos- 
sible now and in foreseen future. The aims of the PNC 
measurements in Yb are different @: (i) to study the 
ratio of the PNC amplitudes for different isotopes, and 
(ii) to measure nuclear spin-dependent PNC effects, such 
as the effect of nuclear anapole moment. The study of 
the PNC for a chain of isotopes does not require atomic 
calculations and can deliver useful information about ei- 
ther neutron distribution or new physics beyond standard 
model (see, e.g. fl4 - fl7j |). The extraction of the anapole 
moment from the measurements does require atomic cal- 
culations, however, high theoretical accuracy is not crit- 
ical here. 

Ytterbium is a very good candidate for both types of 
the experimental studies. It has seven stable isotopes 
with large difference in neutron numbers AN max = 8. 
Two of the isotopes, 171 Yb and 173 Yb, have non-zero nu- 
clear spin provided by valence neutron. This is especially 
interesting since it allows one to measure the strength of 
the neutron-nucleus PNC potential [l2[ (the anapole mo- 
ment has been measured only for the 133 Cs nucleus which 
has valence proton). 

Calculations of the spin-independent PNC in ytter- 



bium were performed in Refs. [H, [HI EH ■ Calculations of 
the spin-dependent PNC were reported in [2(J H3 ■ The 
results of |20j, |21| for the spin-dependent PNC amplitudes 
are presented as tables of reduced matrix elements of the 
spin-dependent weak interaction for different hyperfine 
transitions. This, in our view, leads to some difficulties 
in interpretation. Reduced matrix elements (RME) are 
very convenient for intermediate calculations. However, 
presenting final results in a form of RME may lead to 
confusion due to their unnatural symmetry properties: 



(F a ,a\\H\\F b ,b) = (-1 



>F a -F b 



(F b ,b\\H\\F a ,a) 



(1) 



Here F a is the total momentum of the state a, asterisk 
means complex conjugation which in the case of PNC 
amplitudes means the change of sign. There is an ap- 
parent disa gree ment between the signs of different RME 
in [20| and 21|. The most likely exp lanation for this in 
our view is that the authors of 20] and [21] presented 
different RME, say (a||#||&) in gjj and (b\\H\\a) in gj. 
At least all sign differences follow strictly the rule ([1]). 

Strictly speaking, the sign of an amplitude is not de- 
fined (since a wave function may be multiplied by an 
arbitrary phase factor) , only the ratio of two amplitudes 
between the same states has definite sign. Neither of the 
work [20l |21| provides a link between the spin-dependent 
(SD) PNC amplitudes and spin-independent (SI) PNC 
amplitudes calculated earlier in [18l . Il9j . This means that 
it is hard to say whether the spin-dependent effects in- 
crease or decrease a particular PNC amplitude. In other 
words, the sign of the spin-dependent interaction con- 
stants, such as the anapole moment, cannot be extracted 
from the measurements when using the calculations of 
[2(1 or (21| and no additional assumptions (note that the 



apparent disagreement between the signs of the ampli- 
tudes in [2(| and [2l[ shows that any guesswork about 
the relative signs of the SI and SD amplitudes is unreli- 
able). 

To avoid this problem, in present paper both spin- 
independent and spin-dependent PNC amplitudes are 
calculated simultaneously using the same procedure and 
the same wave functions. In this approach the relative 
sign of the amplitudes is fixed. This allows for unam- 
biguous determination of the sign of the spin-dependent 
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contribution. The constant of the spin-dependent inter- 
action can be expressed via the ratio of the two ampli- 
tudes. This brings an extra advantage of more accurate 
interpretation of the measurements. The accuracy of the 
calculations for the ratio of the PNC amplitudes is higher 
than that for each of the amplitudes. This is because the 
amplitudes are very similar in nature and most of the 
theoretical uncertainty cancels out in the ratio. 



II. THEORY 

Hamiltonian describing parity-nonconserving electron- 
nuclear interaction can be written as a sum of spin- 
independent (SI) and spin-dependent (SD) parts (we use 
atomic units: K = \e\ = m e = 1): 



PNC 



Hi 



si 



(2) 



where Gf ~ 2.2225 x 10 14 a.u. is the Fermi constant 
of the weak interaction, Qw is the nuclear weak charge, 

a and 75 arc the Dirac matrices, J is the mi- 



en 







clear spin, and p(r) is the nuclear density normalized to 
1. The strength of the spin-dependent PNC interaction 
is proportional to the dimensionless constant h which 
is to be found from the measurements. There are three 
major contributions to u arising from (i) electromagnetic 
interaction of atomic electrons with nuclear anapole mo- 
ment, (ii) electron-nucleus spin-dependent weak interac- 
tion, and (iii) combined effect of spin- independent weak 
interaction and magnetic hyperfine interaction (see, e.g. 
[Tl|). In this work we do not distinguish between differ- 
ent contributions to x and present the results in terms 
of total vc which is the sum of all possible contributions. 

Within the standard model the weak nuclear charge 
Qw is given by [H| 



-0.98777V + 0.0716Z. 



(3) 



Here N is the number of neutrons, Z is the number of 
protons. 

The PNC amplitude of an electric dipole transition 
between states of the same parity |i) and |/) is equal to: 



El 



PNC 
fi 



E 



(/|d|n)H#PNc|*> 



Ei — E n 
{f\H?Kc\n){n\d q \i) 



Ef — E„ 



(4) 



where d = — e J2i r, i i s the electric dipole operator, \a) = 
\J a F a M a ) and F = I+J is the total angular momentum. 

Applying the Wigner-Eckart theorem we can express 
the amplitudes via reduced matrix elements 



El 



PNC 
fi 



-M f 



\JfEf\\dpuc\ 



F f 1 Fi 
-M f q Mi 

JiFA. 



(5) 



Detailed expressions for the reduced matrix elements of 
the SI and SD PNC amplitudes can be found e.g. in 
Refs. [22j and [23]. For the SI amplitude we have 



(J / ,^ / ||dsi||Ji,^> = (-l) J+F<+ ^ +1 

Ji Jf 1 



x^(2F / + l)(2F l + l)| -» £ j 
x y-^ \( J fW d W n > J n)(n, J n \\Hsi\\Ji) 



(6) 



Ei — E n 

+ (Jf\\Hsi\\n, J n )(n, J n \\d\\Ji) 
Ef — E n 

= c(J>,J>,F i ,Ji)^. 

Here c(Ff,Jf,Fi,Ji) is the angular coefficient and the 
sum over n, E'^ does not depend on Ff or Ff. 



E' 



E 



(J f \\d\\n, J n )(n, J 71 \\HsjM 



Ei — E n 
(J f\\H S i\\n, J n )(n, J n \\d\\Ji] 



(7) 



Ef — E„ 

For the SD PNC amplitude we have 

G F 

(Jf,Ff\\dsT>\\Ji,Fi) = -j=x 

x + 1)(27 + 1)(2F, + l)(2F f + I)/ 'I 



I I Fi 
(Jf\\d\\n, J n )(n, J n \\ap\\Ji} 



J n Jf 1 
Ff Fi I 



En — Ei 



(8) 



+ (-1) 



Fr-Fi 





Jf 




/ Jn 


Ji 1 1 


{t 


I 


A} 


\Fi 


Ff I ] 



(Jf\\ctp\\n, J„)(n, J n \\d\\Ji) 



E n — Ef 

In the case of the 1 So— ^ 3 Di transition these expres- 
sions can be significantly simplified. Substituting Fi = I, 
Ji = 0, Ff = 1,1 ±1 = F, J f = 1, J n = 1 we have for 
the PNC amplitudes (z-components) Ep i! p, of the tran- 
sitions between specific hfs states of 171 Yb (I = 1/2) 



Ei i 

2 ' 2 ' 



El 3 

2 ' 2 ' 




Similar expressions for 173 Yb (I = 5/2) are 
1 



Ef, 3 

2 ' 2 ' 



Es 5 
2 ' 2 ' 



Es 7 

2 ' 2 ' 




98 



3375 
3l5 

T 

G3 



E"x, 



E"x, 
E"x. 



(9) 
(10) 

(11) 
(12) 

(13) 
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Here E" is the part of the SD PNC amplitude which is 
independent on F t and Ff. 



E" 



72 ^ 



(J f \\d\\n, J n )(n, J n \\ap\\Ji] 

En — Ej 



(Jf\\ocp\\n, J n )(n, J n \\d\\Ji] 



E n - Ej 



(14) 



Note that if at least two PNC amplitudes are measured 
then the value of n can be expressed via the ratio E" /E' 
of the calculated SD and SI PNC amplitudes. This ra- 
tio is much less sensitive to numerical uncertainties than 
each of the amplitudes. The amplitudes are very similar. 
Therefore, a greater part of the numerical uncertainty 
cancels out in the ratio. For example, if amplitudes © 
and (1101) are measured then 



V3 



V2R+1 
S(R - y/2) ' 



(15) 



where R = Ei i z /Ei^ z and S = E"/(E'Q W ). The ra- 
tio of theoretical amplitudes E"/E' is significantly more 
stable in the calculations than each of the amplitudes. 

III. CALCULATIONS 

We consider ytterbium as an atom with two valence 
electrons above closed shells and use the combination of 
the configuration interaction and many-body perturba- 
tion theory (CI+MBPT, [|1]) to perform the calcula- 
tions. The calculations are very similar to our previous 
calculations of ytterbium polarizabilities [25[ . Below we 
briefly describe the procedure emphasizing some minor 
differences. 



CI+MBPT method 



and E is the correlation operator. It represents terms in 
the Hamiltonian arising due to virtual excitations from 
atomic core (see Ref. [13, El| for details). E = cor- 
responds to the standard CI method. Ei is a single- 
electron operator. It represents a correlation interaction 
(core-polarization) of a particular valence electron with 
the atomic core. E2 is a two-electron operator. It repre- 
sents screening of the Coulomb interaction between the 
two valence electrons by the core electrons. We calcu- 
late E in the second order of the MBPT. We use a B- 
spline technique [27J to construct a complete set of single- 
electron orbitals. We use 40 B-splines in a cavity of ra- 
dius R — 40 0,3 and calculate the eigenstates of the V N ~ 2 
DHF Hamiltonian up to the maximum value of the angu- 
lar momentum l max — 5. The same basis is used in com- 
puting E and in constructing the two-electron states for 
the valence electrons. 40 out of 60 lowest-energy states 
for every I up to l m ax = 5 are used to calculate E and 
16 lowest states above the core are used for every / up to 
lmax — 4 to construct the two-electron states. 

The two-electron valence states are found by solving 
the eigenvalue problem, 



H B ** V = E V * V 



(19) 



using the standard CI techniques. Calculated and experi- 
mental energies of a few lowest-energy states of Yb can be 
found in Ref. [25| . The pure ab initio energies are already 
close to the experimental values. However, for improving 
the accuracy further, we re-scale the correlation operator 
Ei by replacing Ei in the effective Hamiltonian (|16[) in 
each partial wave s,Pi/2>J>3/2) • • • by / Ei. The rescal- 
ing factors are f s = 0.875, f p = 1.268, fd = 0.935, and 
ff = l. These values are chosen to fit the experimental 
spectrum of Yb. Some differences in scaling parameters 
compare to what was used in Ref. [25[ is due to the fact 
that in present work we have fitted exactly the energy of 
the state while in j25| we fitted the energy of the 
3 PJ state. 



The effective CI+MBPT Hamiltonian for two valence 
electrons has the form 

H eB = hi(n) + k(r 2 ) + h 2 {n,r 2 ), (16) 

where hi is the single-electron part of the relativistic 
Hamiltonian 

Ze 2 

h-L = cap + (/3 - l)m e c 2 + V N ~ 2 + Ei, (17) 

r 

and h,2 is the two-electron part of the Hamiltonian 

e 2 

h 2 {r u r 2 ) = r + E 2 (n,r 2 ). (18) 

Fi - r 2 | 

In these equations, a and are the conventional Dirac 
matrices, V N ~ 2 is the Dirac-Hartree-Fock (DHF) poten- 
tial of the closed-shell atomic core (N — 2 = 68, Z = 70), 



B. Dalgarno-Lewis and RPA methods 

Matrix elements are found with the random-phase ap- 
proximation (RPA) [H, El 

El ma = (*v\\? + fiV N -*\\* v> ), (20) 

where 5V N ~ 2 is the correction to the core potential due 
to core polarization by an external field /. In present 
calculations / represents either external electric field, SI 
weak interaction or SD PNC interaction. 

Computing PNC requires summing over a complete 
set of two-electron states (see, e.g. Eq. (jU). We use 
the Dalgarno-Lewis method 30] for the summation. In 
this method, a correction 5^ v to the two-electron wave 
function of the state v is introduced and the amplitude 
is reduced to 

A vw = (*tf„||/i| !*„,). (21) 
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TABLE I: Magnetic dipole hyperfine structure constants A 
(MHz) for the 3 Pf and : PJ states of 171 Yb, comparison with 
experiment. 



State 


Calculations 


Experiment" 




4460 


3958 


Ipo 


-819 


-1094 



"Reference 



The correction S^f v is found by solving the system of 
linear inhomogeneous equations 



(H ee -E v )8V v = -(f 2 + 8V 



(22) 



Here /i and f 2 are electric dipole and PNC interaction 
operators (/i = d, f 2 — -ffpNC or vice versa). 



C. Accuracy of the calculations 

Accuracy of very similar calculations of polarizabili- 
ties of ytterbium were studied in detail in our previous 
work 123 and were found to be about 5%. However, 
we cannot claim the same accuracy for present calcu- 
lations due to two important differences. First, there 
is a resonance contribution to the PNC amplitude in- 
volving the 1 P° state. Energy interval between the 3 Di 
and 1 P° states is very small. Its experimental values is 
just 579 cm -1 . The term in (j4|) involving the X P° state 
gives more than 80% of the total PNC amplitude. Even 
very accurate calculations may give significantly differ- 
ent value of small energy interval which would lead to 
large error in the PNC amplitude. One way around this 
problem is to separate the resonance term from the rest 
of the sum and use the experimental energy for the de- 
nominator. We use a technically more simple procedure. 
We have rescaled the correlation operator S to fit the in- 
terval exactly. As a result, the contribution of the error 
in the energy denominator to the error in the amplitude 
is small. 

Another important difference of present calculations 
from those of Ref. [25[ is that we need to calculate matrix 
elements of weak interaction which are sensitive to the 
wave functions on short distances. A way to test the 
wave functions on short distances is to calculate hyperfine 
structure (hfs) constants. 

Calculated and experimental values of the magnetic 
dipole hyperfine structure constants A for the 3 P° and 
1 Pf states of 171 Yb are presented in Table |TJ The first 
calculated hfs constant is larger than the experimental 
one by 13%, the second is smaller by 25%. The rea- 
son for the calculated hyperfine constant of the 1 Pf state 
to differ significantly from the the experimental value is 
the same as for the electric dipole transition amplitude 
between this and ground state - the admixture of the 
4f 13 5d6s 2 (7/2,3/2)? state at E = 28857cm" 1 (see Ref. 
[25( for details). This admixture is small. However, it 



TABLE II: PNC amplitudes (z-components) for the 
|6s 2 , 1 S ,Fi) -> |6s5d, 3 Di,F 2 ) transitions in 171 Yb and 173 Yb 



in units of E'Qw and 10 ieao- 



A I Fi F 2 



PNC amplitude 



units: E'Q\ 



units: 10 ieao 



171 0.5 0.5 0.5 

0.5 1.5 

173 2.5 2.5 1.5 

2.5 2.5 

2.5 3.5 



-(l/3)(l-0.0161x) 
^279(1 + 0.0081x) 



- y4745 (l - O.Olllx) 
- ^5721 (1 - 0.0032x) 
^2/21(1 + 0.0079;*) 



6.15(1 -0.0161^) 
-8.70(1 +0.0081^) 

5.61(1 - 0.0111^) 
9.18(1 - 0.0032x) 
-5.81(1 + 0.0079^) 



can change hfs of the 1 P° state significantly due to the 
large hfs in the admixed state. In contrast, it cannot 
change that much the weak matrix element between the 
1 P° and 4/ 14 5c?6s 3 Di states. This is because the transi- 
tion between the 4/ 13 5d6s 2 and the 4/ 14 5d6s 3 Di states 
in zero approximation is the 6s — > Af transition and cor- 
responding weak matrix element is zero. Therefore, poor 
accuracy for the hfs of the *P° state is not a good indi- 
cator for the accuracy of the PNC calculations. A 13% 
error in the hfs of the 3 P? state gives a more realistic 
estimate for the uncertainty. 

We stress that the uncertainty in the ratio of SD and 
SI PNC amplitudes (E"/E r ) is significantly lower. Tests 
show that this ratio is three to five times less sensitive 
to the variation of the calculation procedure than each 
of the amplitudes. We believe that 10% is a reasonable 
estimate for the theoretical uncertainty for this ratio. 



IV. RESULTS 

Calculations give the following value of the spin- 
independent PNC amplitude of the 1 So — > 3 T>i transition 
in ytterbium: 



£ SI-PNC = L12 g x 10 -11 Q 



(23) 



This corresponds to the following value of the reduced 
matrix element 



E' = 1.945 x 10" 



iea,Q. 



(24) 



The electron (^-independent) part of the reduced matrix 
element of the spin-dependent PNC amplitude is found 
to be 

E" = 3.648 x 10" n iea , (25) 

The effect of different nuclear size for 171 Yb and 173 Yb 
is only 0.1% for both SI and SD PNC amplitudes. It 
is neglected in (gSJ), (EH) and d^HJ) • We use Fermi- type 
distribution for nuclear density p with nuclear radius 
R N = 6.35 fm for 171 Yb and R N = 6.37 fm for 173 Yb [13]. 
The ratios of the SD and SI PNC amplitudes are 



E"/(E'Q W ) 
E"/{E'Q W ) 



-0.0198(20) for m Yb, 
-0.0194(20) for 173 Yb. 



(26) 
(27) 
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The difference in these values is due to different weak 
nuclear charge Qw {Qw — —94.75 for 171 Yb and Qw = 
—96.72 for 173 Yb). The difference is within numerical 
uncertainty. 

The results for the specific PNC amplitudes between 
different hyperfine structure states of m Yb and 173 Yb 
arc presented in Table HU These results are obtained by 
substituting (J24J) and ([25]) into ©, (10]), (JTU), (12]) and 
([T3"]) . The numerical factors before x are proportional to 
(E" /E'). The theoretical uncertainty for these factors is 
about 10%. The expressions from the table or equations 
(0), (HO]), (HU), (US), (HU) together with (24]) and (25]) can 
be used to extract the value of x from the measurements. 
For example, for m Yb Eq. (15]) becomes 

^=88(9) 1 + ^ R . (28) 



A. Comparison with other calculations 

Calculations of the spin-independent part of the 
PNC amplitude for ytterbium were performed before in 
Refs. P, IH, [lil- The spin-dependent amplitudes were 
calculated before in Refs. [2(J Hl[. ^ is convenient to 
compare the results in terms of E' (7]) and E" (T4"]) 
since these values are the same for all hfs transitions. 
Refs. pj [H, [l9| present the values of the z-component 
of the SI PNC amplitude. Refs. [H EJ present reduced 
matrix elements of the SD PNC interaction for each hfs 
transition. Corresponding values of E' and E" can be 
easily extracted from this data using formulas of present 
paper. The absolute values of the amplitudes are pre- 
sented in Table IIIII We have excellent agreement for E' 
with DeMille [l[ and Porsev et al [l8[ while the result of 
Das [ill is about 30% smaller. We have good agreement 
with both Porsev et al [2l| and Singh and Das [2(| for 
E". The difference between Ref. ]W and our result is 
12% which is within our uncertainty. The difference be- 
tween our results for E" and those of Ref. t 2(| is even 
smaller. But this is probably accidental. Note however 
that we agree with Singh and Das [2(| on small, practi- 
cally negligible change of E" from m Yb to 173 Yb while 
Porsev at al [2l| report a 3% increase. Such increase has 
no physical explanation and must be a numerical effect. 
In our experience such effect can be a result of just one 
RPA iteration after a change of nuclear radius from 171 Yb 
to 173 Yb. Further iterations kill the difference. However, 
it is up to the authors of [21| to explain their results. 

Note again that in Table lllll we present only the abso- 
lute values of the amplitudes, ignoring their signs. This 
is because the sign of an amplitude is not fixed and has 
no physical meaning. However, the relative sign of the SI 
and SD PNC amplitudes is not arbitrary. The SD depen- 
dent part of the PNC amplitude must either increase or 
decrease the transition amplitude depending on the sign 
of x. It is important to know the relative sign of the am- 
plitudes to be able to extract the sign of x from the mea- 



TABLE III: Spin-independent (E 1 ) and spin-dependent (E") 
parts of the PNC amplitude (reduced matrix elements) for 
the |6s 2 , 1 So) — ► |6s5d, 3 Di) transition in ytterbium, compari- 
son with other calculations. The signs of the amplitudes are 
omitted. 





E 1 


E" 


(10- n iea ) 




10 _11 iea 


171 Yb 


173 yb 


DeMille [1] 


1.9 






Porsev et al 


1.97" 


4.13 6 


4.27 6 


Singh and Das 


1.33 c 


3.68 d 


3.67 d 


This work 


1.95 


3.65 


3.64 



"Ref. [3, 6 Ref. [2lJ, c Ref. [TjJ, d Ref. 



surements. To fix the relative sign of the two PNC am- 
plitudes one should calculate them using the same wave 
functions. This is how it is done in present work (see Ta- 
ble M and formulas ©, (ITJ, (TTJ, (TJJ), ([15])). Another 
important advantage of the simultaneous calculation of 
both amplitudes is that x can be expressed via the ratio 
of the amplitudes. This ratio has much smaller theoreti- 
cal uncertainty than each of the amplitudes (see previous 
section for discussion). 

Unfortunately, both previous calculations of the SD 
PNC amplitude in Yb [2^, [2l| do not compare their re- 
sults with the earlier calculations of the SI PNC ampli- 
tudes [H, EH which could be performed using different 
wave functions. This leads to the uncertainty of the rel- 
ative signs and larger errors in the ratios of the SI and 
SD amplitudes which are needed to extract the value of 
x from the measurements. 

B. Comparison with experiment 

The result of the measurement of the PNC amplitude 
of the 1 So — > 3 Di transition in 174 Yb reported in [6| reads 

|£ PN C| = 8J ( L4 ) x i - 10 ea . (29) 

Assuming a 13% theoretical uncertainty and substituting 
weak nuclear charge Qw = —97.71 we get from (j2"3"j) the 
following theoretical value for the amplitude 

|£ PNC | = 11.0(1.4) x 10- 10 ea . (30) 

The values of (|29p and ([317]) agree within the declared 
uncertainty. 

To measure the constant of spin-dependent PNC inter- 
action [x) more accurate measurements are needed for 
171 Yb or 173 Yb. The work is in progress at Berkeley @. 



V. CONCLUSION 

We present simultaneous calculation of the spin- 
independent and spin-dependent PNC amplitudes of the 
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6s 2 1 So — > 6s5d 3 Di transition in ytterbium. The re- 
sults are to be used for accurate interpretation of fu- 
ture measurements in terms of the parameter of the spin- 
dependent PNC interaction x. Both, sign and value of 
k can be determined. Theoretical uncertainty is at the 
level of 10%. 
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